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ABSTRACT 

Propagation of waves in curved duets and pipes belong to the class of 
motion 'which is characterized by wave patterns totally different from 
those known in straight duets or in unlimited space. The curvillinear 
boundaries are responsible for the appearance of a continuous standing 
radial wave which in turn affects the transmitted tangential waves . 

The purpose of this paper is to solve the problem of propagation of long 
acoustic waves in slightly and sharply bend duets. This problem has 
been only partially analyzed by various authors . In this study two 
acoustic systems are considered which allow determination of the basic 
modes of motion and describe the transition and distortion of plane 
waves as they propagate down the curved channel . A detailed study of 
the behavior of waves in junctions between straight and curved ducts is 
also given. Solutions and expressions for principal modes of the wave are 
obtained by using the linearized equation of motion solved for its 
characteristic values. This original approach required a novel use of 
Bessel functions to determine the characteristic values of the steady 
and the decaying fields of motion. 


INTRODUCTION 

#• 

The propagation in curvilinear waveguides or bent ducts is difficult 
to analyze, and mathematical models developed to date are complicated. 
Perhaps that is one of the reasons why relatively few papers are avail- 
able on this subject. 

The first recognition of the problem of propagation of pressure waves 
in a curved conduit, as distinct from the motion in a straight line, 
was formulated by Rayleigh (Ref. l) , In a short, brilliant expose he 
demonstrated that long waves in a curved pipe of infinitesimal cross 
section behave exactly as in a straight pipe. The curvature of the 
pipe is of no importance. His analysis is based on the linearized 
equation of motion and conclusions are not valid if the cross section 
of the pipe is finite. 
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Subsequent research, of wave motion in curvilinear ducts is almost 
exclusively analytical and the majority of papers treat the propaga- 
tion of electromagnetic waves in curved ducts . Only a few discuss the 
propagation of sound waves . Interestingly enough, along with analyti- 
cal formulation of the behavior of waves in bends there appeared a 
series of papers dealing with mathematics needed to solve the physical 
problem of waveguides . 'This parallel effort indicates that the solution 
of the problem of bent waveguides requires mathematical formulations and 
techniques not generally available, A number of papers published between 
1939 ana 1969 did contribute to the problem, Buehholz (Ref, 2), using 
separation of variables obtains a solution for propagation of electro- 
magnetic waves in slightly bent waveguides of infinite length. He cal- 
culates an. expression for a wave propagation constant and draws several 
general conclusions about behavior of waves in bends. Buehholz 's paper 
is the first contribution to the problem of progressive waves. The 
problem of the infinitely long bend was also treated by Krasnushkin 
Vhex . 3) , Grigor'yan (Ref, 4), and others. Krasnushkin approaches the 
problem by the method of separation of variables but in view of mathe— 
maoicai difficulties proposes a perturbation method and treats the sim- 
plilied ease of the slightly bent tubes, Grigor'yan solves the differen- 
tial equation by expansion of the cross product of Bessel functions into 
a Taylor series. He tries to obtain an algorithm of sufficient generality 
uC/ be applicable to all possible impedances of the waveguide walls . 

The method is only partially successful. Grigor'yan obtains correct 
genej.al information on amplitude and distribution of the radial os- 
C ^ ilat ^ ons basic mode wave number does not verify the differen- 

c.i_al equation except for the Raleigh case of very narrow pipe. 

Among the few treatises on propagation of sound in bends are 
those by ^ Miles (Refs, 5 and 6) in which he establishes an analogy between 
propagation oi sound and an electrical transmission line are most im- 
portant. The method is then applied to a right angle joint of rec- 
tangular tubes. 


The work of Miles was checked by Lippert (Refs. J to 9). Lippert 
presents an experimental study of sound wave propagation in mitered 
uends of various angles . The experiments were conducted over a wide 
- ange of frequencies and show that long waves in mitered bends propa- 
gate with insignificant reflections. 


. present analysis considers propagation of long acoustic waves in 
a circular bend after Rostaf inski (Ref, 10), Established will be the radial 
' 'f' e a -8 u ' ! - ax distribution of the 'vibrational velocities. The mathemat- 
j rea '" m ‘- n f of '-he problem utilizes the method of separation of variables, 
solutions and expressions for principal modes of the wave are obtained by 
J~xng one linearized equation of motion solved for its characteristic values. 
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GENERAL SOLUTION 


Determination of Eigenvalues 


The linearized wave equation in cylindrical coordinates is known to 
he separable, and since the superposition of solutions Is allowed the 
general solution may he written 


i (tot+a ) f 
e (. a 
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v v 
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v v 


(kr )[] dv 
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where C is a set of points in the complex plane , to he determined in 
order to satisfy the boundary conditions. 

The solution for v = 0 might have been included in the integral 
terms but in order to show more explicitly the linear dependence on 9 
of this solution, it has been written separately. 


To satisfy the partial differential equation and the boundary condi- 
tions for perfectly rigid circular walls, a characteristic equation will 
be found whose roots will be the characteristic values of the problem: 
a set of values of the separation parameter v which will yield a non- 
trivial solution of the problem. Differentiating equation (l) with 
respect to r and equating to zero we obtain for the two circular 
boundaries R^ and R^ 


v # 0 (2a) 


v = 0 (2b) 


<< 1, with the steep slope 
of Jq(kr) over the range 
range under consideration. 

Consequently, the solution v = 0 cannot be considered. To 
evaluate the v m 's we expand the J\J, and Jl v in terms of increas- 
ing powers of the argument (akRq) = (kR2) and (kRq) . If we limit the 
expansion in the first approximation to the first term we obtain 
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In the present case where 0 < (kr) < (kR^) 
of Y^(kr) and relatively moderate increase 
(kRq) to (kRg), there is no solution in the 
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Solution v - 0 has teen already rejected. Therefore, the only 
acceptable solution must satisfy the equation 

v-l -v-l 
a = a 

which may be put in the form 
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Beuter approximations will be given by the second and following terms of 
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expansion of J’ and J» The result is 
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Tliere is an infinite set of pure imaginary roots v given by equa- 
tion ( 5 ) and one single real root.. The obtained eigenvalues allow us 
to write a complete solution in which only the integration coefficients 
in m are to be determined by proper boundary condition at the inlet 
and outlet of the duct. This determination is a relatively simple matter. 


TYPICAL PHYSICAL SYSTEMS 


Infinite Bend 


Suppose an infinitely long circular duct (Fig. la) for which the 
closest physical example is a tightly wound coil of which the pitch is 
negligible compared to the radius of curvature of the duct . The far 
end boundary condition for an infinite duct is that no reflection of 
waves must be considered. 

The integrated equation for the velocity potential <j>(r,0,t) for an 
infinite bend may be written directly by substitution of the derived ex- 
pressions for the integration constants 
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The propagation of waves' in a curved duct is profoundly influenced 
by the curvature of the' conduit. The amplitude of vibrational velocities 
is a function of "both inner radius R-j. and of the radii ratio a = R 2 / R l* 
The tangential vibrational velocities are almost exactly inversely pro- X 
portional to radius, that is, velocities follow the distribution of a 
potential vortex. The radial vibrational velocities are approximately 
two orders of magnitude smaller than the tangential velocities. 


The phase velocity of waves is # also affected by bending of ducts, 
i'or a bend, the phase velocity is 0 = w/v , as obtained from 

(wt - v o 0), while for a straight duct it Is x = to/k, as given by 

(ait - kx). To compare _ the two velocities we average the tangential 

phase velocity, 0r = s over the duct width and obtain Is! = ai/(v^)R m 

The ratio of the two velocities is 0 mean 
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Clearly , the phase velocity in bends is always 'higher than in straight ducts 


The attenuated tangential vibrations which characterize change from 
motion of plane waves in straight duct to motion in a curved duct are 
examined in the next two figures . Figure 2 gives results of a sample 
calculation illustrating the behavior of those oscillations for a duct 
of radii ratio a = 2. The vibrations are basically of low amplitude. 
Even^ close to the piston, at 0 = tt/1 6, they are one order of magnitude 
smaller than vibrations of the nondamped, propagating wave. The radial 
distribution of these oscillations changes significantly with wave angular 
position in the duct. At 0 = ir/k these oscillations are reduced to 
a very low level and are nearly uniform across the width of the duct, 
tiguie 3 shows the same oscillations calculated for three different duct 
widths but with a single angular position of 0 = ir/k. The curve for 
a - 2 was taken directly from figure 2 for comparison purposes . The 
two other curves indicate that the decaying oscillations are much more 
pi enounced and extend farther when Induced In wider ducts. 


The radial vibrational velocities, which characterize motion of 
waves in curved ducts, for long waves, are of low amplitude. The per- 
manent, standing oscillations are shown in figure 4. They are calculated 
i oi^ duct radii ratios of 2, 3, and 4. Generally the amplitude of these 
us ex nations is low, approximately two orders of magnitude smaller than 
cne tangential velocities. The radial distribution is characterized by 
trie lack of symmetry. The maxima of curves are shifted toward bend's 
inner wall. 

This phenomenon is even more pronounced in the case of the nonpro- 
pagating, damped, radial oscillations at the curved duct Inlet. The 
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amplitude and radial distribution of these oscillations is shown on 
figure 5 for three values of parameter a and 8 = n/k. For a = 2 the 
permanent and the vanishing oscillations are of approximately the same 
amplitude. For a = 3 and k the decaying oscillations are about twice 
as large as the' radial oscillation of the propagating wave. Therefore, 
the process of decay is much slower when ducts are wide. 


Radians Bend Followed by a Straight Infinite Duct 

To analyze the motion in the system bend-straight duct, it is neces- 
sary to solve 


1 a tf) 8(j) 


3<j> _ 9tf> 


and 

r 38 3x 3r 3y 


at 


0 = 8 , 


x = 0 


that is, the continuity conditions at the junction bend-straight tube, 
of a system illustrated on figure lb. 

Solving the differential equations and evaluating the constants of 
integration we obtain the velocity potential 
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where the term e 0 2 results from evaluation of a set of integration 
constants at x = 0 which is at 0 = 8g and not at 0=0 and could 
be avoided should the coordinate 0 of the bend be counted clockwise 
from the junction bend-straight duct. 8 m is a set of coefficients 
resulting from matching of vibrational velocities at the junction. 
are smaller than one . 

The derived solution is very general. It applies to radii ratios of 
the bend in the range from 1 to 10; it has no limitations as to the 



length, of the bend which may range from a small fraction of it to an 
infinite coil. This wide range of applicability was obtained by retain- 
ing in the solution, terms' which are two orders of magnitude smaller 
than the principal terms. The first two terms of equation (8) are those 
of equation (7) of the wave motion in an infinite bend. The first term 
is a simple wave, the second the attenuated vibrations near the inlet 
section of the bend. The third term represents attenuated vibrations 
j.,n the negative direction of 0, an adjustment of wave motion to discon- 
tinuity at the junction bend-straight duct. This term is identically 
zero at .0 = 0 and increases to a maximum at 6 = 0 2 . 

Ihe velocity potential in the straight, infinite duct is approxi- 
mately 
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where A mn is a set of integrals in m and n. 

The calculation of the vibrational velocities in the bend and in the 
straight duct that follows the 90 degree bend is straightforward. Re- 
sults are reported on figure 6 . 

tIle straight . duct the waves are straightening out relatively fast. 
' 'y X/f '^ ~ 1 there is already an almost straight wave similar to that 
i*hich was generated at the bend inlet. Essentially, the process of 
straightening out of the wave is confined to the straight duct. 

The verification of the simultaneous solution of the equations of 
motion for the bend and for the straight duct is shown on figure 7. 
ihe calculated tangential vibrational velocities for r = Rt and 
r -^2 _ for positions starting at the vibrating piston, through the 
bend, and in the straight duct are taken from figures 11 and 12. The 
rapid changes in velocities by the bend exit are well illustrated. In 
spite of this, the values calculated for the bend match very well the 
values calculated for the straight duct. 
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CONCLUDING REMARKS 

Propagation of long wayes in a two-dimensional system has been 
analyzed. The acoustic approximation has "been used and the He lmh ol z 
equation was integrated for two sets of "boundary conditions. The two 
physical systems taken into consideration are: an infinite "bend ap- 

proximating a coil and a 90 ° "bend followed "by a straight, infinite 
duet approximating a typical industrial piping system. 

The results of the analysis indicate that "bending of a straight 
duct profoundly modifies the propagation of waves in that duct. The 
"bend engenders the following phenomena: 

1. a set of attenuated axial and radial waves which modify the 

plane wave generaged at the duct’s inlet 

2. a continuous radial, standing wave whose radial vibrations are 

sustained by the curvature of the bend. In a straight duct 
these vibrations would be quickly attenuated 

3. a vortex- type distribution of the tangential vibrational veloc- 

ities . 

The presented analysis Is not directly applicable to real flows 
because it is based on a linearized equation valid for acoustical waves 
In stationary medium only. It gives, however, an idealized picture of 
wave's behavior In bends, in general. 
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Figure 3. - Attenuated tangential vibrational velocities for three bends at 0 = tt/4. 



Figure 4. - Standing radial vibrational velocities 
in bends for three R, and three values of a. 
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(a) Tangential vibrational velocity distribution near inlet of an infinite coil. 
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(b) Distribution of tangential vibrational velocities at bend outlet. 
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(c) Straight duct inlet. Distribution of axial vibrational velocities. 
Figure 6. - 90° Bend followed by a straight duct. Distribution of vibrational velocities. 
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Fig ure 7. - Propagation in bend-straight duct system. Tangential vibrational 
velocities at bend's curved walls, at r = Rj and r=R 2 . 
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